Forward-backward tt asymmetry from anomalous stop pair production 
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We analyse a simple Standard Model (SM) extension with only two new light fields: a scalar 
partner of the top i (with mass above mt) and a light neutral fermion x° (with mass of a few GeV), 
coupled to SM quarks via a Yukawa interaction. We show that such model can lead to a significant 
enhancement of the forward-backward asymmetry in tt production at Tevatron via the additional 
it pairs produced from tt' decays. The model satisfies existing constraints on new-physics searches 
both at low and high energies, and could even address the cosmological dark-matter abundance. 
The implications for future searches at the LHC are briefly outlined. 



I. INTRODUCTION 

There are recent experimental indications of an anoma- 
lously large forward-backward asymmetry (FBA) in top- 
antitop pair production at the Tevatron. The asymmetry 
is defined as 

A ti _ cr Ay>0 - OAy<0 s 
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where Ay — y t — yi is the difference in rapidity, and pos- 
itive rapidity is measured in the direction of the colliding 
p. The most significant measurement of A l p B is due to the 
CDF collaboration which after subtracting backgrounds 
and performing the unfolding procedure to recover the 
asymmetry at the partonic level reports [1] 

Af B =0.158 ± 0.074, (2) 

to be compared with the NLO QCD prediction 
(Af B ) S M = 0.058 ± 0.009 [2]. The observed FBA is 
even larger at high invariant masses of the tt system, in 
particular 

A h = Af B{mt{>450GcV) = 0.475 ±0.114, (3) 

while QCD predicts 0.088 ± 0.013 [2]. 

At the same time, the total inclusive top pair produc- 
tion cross-section is currently measured to be [3] 

cr" = (7.50 ± 0.48) pb , (4) 

that is consistent with the most recent theoretical SM 
predictions for this observable: (7.2 ± 0.4) pb [4] and 
(6.4 ± 0.4) pb [5]. Another important experimental con- 
straint is the m t f distribution of the production cross- 
section. As pointed out in [6], the most significant infor- 
mation at high m t i is the one derived from the next-to- 
highest measured bin [7] 

a h = C(7 00 GcV<m tt -<800 GoV) = ($0 ± 37) fb , (5) 
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to be compared with {<Jh)sM = 80 ± 8 fb [5]. 

Several authors have analysed the possibility that the 
large FBA reported by CDF is obtained by an anoma- 
lous ti production mechanism interfering with the corre- 
sponding SM process [6, 8]. The interference of the two 
amplitudes maximizes the possible impact on A l p B while 
minimizing deviations from the SM in a u . However, we 
observe that the present cross-section measurements still 
leave some room for an additional (incoherent) produc- 
tion of ti pairs. This motivates us to analyse a different 
mechanism to enhance A l p B , namely an anomalous pro- 
duction of tt + invisible particles, which do not interfere 
with any SM process. 

At the level of the inclusive measurement, the non- 
standard production (passing the experimental selection 
cuts) can still contribute of up to 13% compared to the 
SM (using the higher SM prediction in [4] as a conserva- 
tive reference value). The percentage of the new contri- 
bution could even rise at high m t f, provided it does not 
exceed 50% of the SM cross-section in ah- Ideally, if the 
anomalous production would carry a 100% FBA, it could 
perfectly accommodate the measured value of the total 
FBA without violating the bound on the cross section. 

The production of the ti + invisible final state can 
be obtained by the pair production of "top partners" - 
particles decaying into a top quark and light invisible 
states- that are naturally expected in several SM exten- 
sions. If the mass difference between the top partner and 
the top is sufficiently small, the missing energy carried by 
the invisible states is small and the ti pairs thus produced 
would pass the experimental cuts applied to identify ti 
pairs in the SM. A prototype of such scenario would be 
the fourth generation up quark or a vector-like fermionic 
top partner, decaying into a top and one or several light 
invisible particles. However, colored fermions have large 
QCD cross-sections with a vanishing FBA contribution at 
leading order: this makes them them unattractive can- 
didates for our purpose. On the contrary, a scalar top 
partner (t) of mass around 200 GeV decaying into a top 
and a single invisible particle (x°) is still perfectly al- 
lowed [9]. This is because the QCD production cross- 
section for scalar particles proceeds mainly via p-wave 
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and thus vanishes at threshold [10]. 

The stop and a light neutralino are naturally present 
in supersymmctric extensions of the SM, making this 
scenario particularly attractive. As we will discuss be- 
low, in the minimal supersymmetric extension of the SM 
(MSSM) the FBA generated by this mechanism is van- 
ishingly small, but it could become sizable in more gen- 
eral frameworks. Without assuming a specific model, 
here we adopt a phenomenological bottom-up approach: 
we assume t and x° 1° be the only light relevant non- 
standard particles, and we determine the nature of their 
interactions mainly by looking at phenomenological con- 
straints. As a remnant of i?-parity in the MSSM, we also 
assume a discrete Zi symmetry under which only these 
non-standard states are charged, such that they can only 
be produced in pairs and such that x° is a stable particle. 



II. IDENTIFICATION OF THE MODEL 

In order to generate a large A*p B , the differential par- 
tonic cross-section u(pi)u(p 2 ) — > t(j>i)P (p'2) should ex- 
hibit a large t-odd dependence, 1 where t = (p[ — pi) 2 ■ 
This dependence is not generated by the leading QCD 
contribution to uu — > tft . 

A first approach to generate a sizable t-odd dependence 
m uu — » is via heavy mediators (in either J or t chan- 
nel). Integrating out the heavy mediators leads to an 
appropriate set of higher-dimensional effective operators 
coupling the up quarks to the tP pair. Up to canonical 
dimension six, the relevant operators are 

unfit, U'yzuPi, vr/nuvd^t, uj^uPd^f . (6) 

In this case the required t-odd dependence can only ap- 
pear in the numerator of the cross section, but the kine- 
matical structure of all the above operators is such that 
they do not generate it. A non-vanishing t dependence 
can appear only introducing operators of dimension 7 
or higher, which are expected to be more suppressed by 
naive dimensional analysis. Since the contributions with 
a t-odd dependence should dominate the cross section, 
which receives also a sizable contribution from ordinary 
QCD, this approach appears to be highly contrived. 

A second approach is to assume light t-channel medi- 
ators for uu — > W. In this way the t-odd dependence 
is naturally induced by the t-channel propagators and 
gets more pronounced the lighter the mediators. The 
lightness of the mediators is also required to generate a 
sizable cross-section (larger than the QCD induced one) 
while maintaining perturbativity of the associated media- 
tor couplings to (light) quarks and t. The mediators need 



The uu initial partonic state is chosen because it largely domi- 
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to be neutral and, more generally, SU(2)l x U(l)y sin- 
glets, to avoid the bounds on new light states from LEP. 
Assuming they couple trilinearly to u and t, they should 
also be fermions. Since we have already introduced a 
light neutral fermion in order to account for t — > tx , the 
simplest possibility is to assume x° itself to be the me- 
diator. We are thus led to consider the following simple 
Lagrangian 

C = C SM + (D^(DH) ~ mlPt + x°(i7uD»)x° 

-rn x X°X° + E (Y q q R t X ° + h.c.) , (7) 

q— u,c,t 

where we have introduced effective fetx couplings for 
all three generations of right-handed up quarks. Having 
assumed x° to be an SU(2)l xU(1)y singlet, this implies 
that t has the same SU(2) lxU(1)y quantum numbers of 
right-handed up quarks. As we will discuss below, beside 
the minimality of the new states introduced, the choice 
of a pure right-handed coupling minimizes the impact 
on electroweak observables and flavour-changing neutral- 
current (FCNC) transitions. 

Following the principle of introducing the minimal set 
of relevant new states, we have assumed x° to have only 
a Majorana mass term. This hypothesis has essentially 
no impact on collider phenomenology (provided m x is 
sufficiently small in order to suppress the production of 
same-sign tops) , while it may be an important ingredient 
if we require x° to be a viable dark-matter candidate (see 
below) . 

With this choice t can be identified with the right- 
handed stop of the MSSM or, more generally, with any 
combination of right-handed up-type squarks. In princi- 
ple, x° could be identified with the bino of the MSSM. 
However, bino couplings are completely determined by 
electroweak gauge symmetries and ultimately turn out 
to be too small to significantly affect the FBA. On the 
other hand, we see no obstacles to accomodate the re- 
quired framework in extensions of the MSSM with extra 
SU{2)l x U(1)y singlets. An important condition in 
order to implement this mechanism in supersymmetric 
extensions of the SM is a sufficeintly heavy gluino mass, 
such that gluino-mediated amplitudes do not provide a 
significant enhancement of the tP production cross sec- 
tion. 

As anticipated, we keep our discussion general 
analysing the phenomenological implications of (7) with- 
out assuming extra model-dependent conditions. In or- 
der to avoid the production of x° pairs in charm decays 
we require m x > itid/2, while for the t state we require 
fiij > m t and we impose the contraints derived in [9] for 
top-partner pair production at Tevatron. 



III. COLLIDER PHENOMENOLOGY 

As an illustration of the resulting phenomenology at 
the Tevatron, we choose = 200 GeV and m x o = 



3 




195 200 205 210 215 
m ; [GeV] 



Figure 1. Tevatron constraints in the m~Y u plane. The in- 
clusive A l p B and a are reproduced within 1 a in the central 
green band. The region below the continuous (dashed) blue 
line is excluded by A h at 95% C.L. (90% C.L.). The region 
above the continuous (dotted) red line is excluded by ah at 
95% C.L. (90% C.L.). 



2 GeV. We simulate our signal pp — >■ W — > ttx°x° & t the 
partonic level with MadGraph/MadEvent 4.4.57 [11] and 
using CTEQ6L1 set of PDFs [12]. For Y u ~ 0.85 we find 
that the total tt' production cross-section reaches 12% of 
the SM ti cross-section, exhibiting a 50% FBA. Assum- 
ing t decays dominantly into t\° (or assuming Y t 3> Yu,c) 
and applying the ti reconstruction cuts, this leads to a 
predicted A l p B within one standard deviation from the 
experimental data in (2). At the same time, with this 
parameter choice the value of the total cross-section is 
within one standard deviation from the measured one, 
using the prediction in [4] as a conservative SM normal- 
ization. Some tension does develop only when comparing 
the high invariant-mass data, where the predicted cross- 
section and FBA tend to be higher and smaller than data, 
respectively. In particular, the predicted <jh, exhibits a 
tension of 1.6 a when compared with (5), while he SM 
tension of 3.4 a in Ah is only relieved to around 1.9 a. 

In general, as long as x° is much lighter than t and 
Y t iu,o the results do not depend on the precise val- 
ues of m x o and Y t . For this reason, we have performed 
a more accurate evaluation of the various constraints on 
the parameter space of the model at fixed m x o and Y t . 
In particular, we have set m x o = 2 GeV and Y t — 4 (the 
largest value to ensure a perturbative behavior in Y t ) and 
varied the remaining relevant parameters, namely Y u and 
nif. Since the condition Y t S> Y u is not always fulfilled, 
we have taken into account the actual t — > t\° branching 
ratio when computing our signals. The resulting com- 
parison with the relevant constraints coming from ti pro- 
duction at the Tevatron is displayed in Fig. 1. 



As shown in Fig. 1, the requirements of a large FBA 
and a small impact on the cross section are partly in 
conflict, especially at high invariant masses. However, 
there is a region of the parameter space where all con- 
straints are satisfied at the 90% C.L. This happens for 
200 < m t ~ < 205 GeV and Y u > 1.5 (part of the green 
area where the dotted line is above the dashed one) . The 
reason for that lies in the interplay of the Y u dependence 
in the x°-mediatcd cross-section generating a large FBA 
and in Br(t — > t\)- In particular, for |V„| < \Y t \ the 
i-even QCD contribution is suppressed by Br{t — > t\) 
more than the £-odd contribution, which enables us to 
obtain a larger FBA value for a given contribution to the 
cross-section. 

The two additional invisible particles in the final state 
result in extra missing energy compared to the SM tt 
production. This might result in (1) observable modifi- 
cations of the various kinematical distributions used in 
the experimental analyses of the FBA (as well as the 
ti production cross-section and the top quark mass) to 
discriminate signal from background; (2) tensions be- 
tween the leptonic, semileptonic and all-hadronic modes 
in the precise measurements of the ti cross-section and 
the top quark mass. Unfortunately, most present top 
quark analyses employ multivariate techniques which are 
difficult to reproduce without the full detector simula- 
tion. Nonetheless, we have estimated the size of such 
effects by comparing the leptons and jets px distribu- 
tions as well as distributions of missing transverse en- 
ergy (E™ lss ) and the scalar sum of transverse particle 
energies (Ht = Yli ^t) m the process pp — > tix°X° ( at 
preferred values of the model parameters) with those of 
the SM ti production. We have considered these distribu- 
tions separately for the leptonic mode with kinematical 
cuts used in [13], the semileptonic reconstruction chan- 
nel described in [14] as well as the all-hadronic mode 
used recently for the top mass measurement in [15]. We 
have simulated both signals using MadGraph, interfaced 
with Pythia 6.4.14 [22] for showering and hadronization, 
and the CDF PGS for detector simulation with a cone 
(dR — 0.4) jet reconstruction algorithm. In the lep- 
tonic mode, we first observe that the fraction of the 
itx X signal passing the kinematical cuts is compara- 
ble to the SM tt signal (the same is true also for the 
semileptonic mode). Regarding the kinematical distri- 
butions outlined above, the effects of the extra x°'s can 
be readily described as a shift in the distributions by 
few GeV. Since such shifts only affect the O(10%) new 
contribution to the total ti sample, the resulting effects 
turn out to be much smaller than the reported system- 
atic and statistical uncertainties in these distributions. 
As a consequence, in di-leptonic semple we expect the 
modifications in the total measured ti cross-section or 
the top-quark mass to be smaller than the present sta- 
tistical uncertainties. In the semileptonic case, the shift 
in the distributions is a bit more pronounced, of the or- 
der of 10 GeV, but the resulting change in the combined 
itx X + it signal sample distributions is still smaller 
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than the associated statistical uncertainties. We have 
also compared the reconstructed transverse leptonic W 

masses (m^ = 2Elf lss E^l — cos(0£ - </> B mi 5 ,)], where 

<pi are the corresponding azimuthal angles) and found 
good agreement between the distributions of the two sig- 
nals. Most important, the distribution peak in the tix°X° 
sample does not appear shifted with respect to the SM 
tt case. Therefore, we expect that resulting effects on 
tt cross-section or top-quark mass measurements using 
the semileptonic mode to be within the present statis- 
tical and systematic uncertainties. Finally, for the all- 
hadronic mode, we find that the requirement of no sig- 
nificant _E™ 1SS reduces the fraction of ttx°X° signal events 
passing the cuts to one third, compared to the SM tt sam- 
ple (assuming an uncertainty in E™ las reconstruction of 
cr B », B3 ~ 0.5yi?rGeV 1/2 [16]). The resulting discrep- 
ancy in the measured total tt production cross-section 
compared to the semileptonic and leptonic modes for pre- 
ferred values of our model parameters would be of the 
order of 7%, that is less than the present uncertainties of 
the individual analyses. 

Recently top-pair production production has also been 
measured at the LHC. The value reported in [17] for the 
total tt cross section in pp collisions at y/s — 7 TeV 
center-of-mass energy is 



"LHC 



(194 ±78) pb, 



(8) 



in agreement with the SM prediction of 158 ± 24 pb [18]. 
For our preferred choice of parameters the corresponding 
tv production cross-section is only around 10 pb [10] 2 , 
still below the present sensitivity of the LHC experi- 
ments. 

The ffi pairs decaying to two up quarks and two x°'s 
contribute to the jets plus missing transverse energy 
(E™ lss ) signatures, searched for both at the Tevatron [19] 
and at the LHC [20, 21]. At present the most sensitive 
search is the ATLAS analysis of two jets plus E™ iss [21], 
based on 35 pb - of data, and in particular the scenario 
A of [21], which is optimized for low-mass squark pair 
production. Here two jets with pt > 140, 40 GeV re- 
spectively, and E™ lss > 100 GeV are required. Addi- 
tional kinematic cuts are also imposed, among which the 
most relevant are m e // = 

Ei=i \pt\ + e t 1ss > 500 GeV 

and E™ lss /m e ff > 0.3. Under these conditions 87 events 
are found passing all the cuts, with a SM background 
estimation of 118 ± 42. Again we have simulated our sig- 
nal contribution using MadGraph, interfaced with Pythia 
6.4.14 [22] for showering and hadronization and the AT- 
LAS PGS for detector simulation with a kr jet recon- 
struction algorithm. For the most interesting region of 
the model's parameter space, we find that the signal 



cross-section passing the pr and E^ lss cuts alone is at 
the level of 1 pb or less, still below the present experi- 
mental sensitivity. 

In principle, also the Tevatron analyses for single top 
production could be used to set bounds on the parame- 
ter space of our model. However, in this case it is very 
difficult to assess the sensitivity. Present experimental 
analyses employ sophisticated multivariate discriminants 
including tight constraints on the missing transverse en- 
ergy distribution, which is required to agree with expec- 
tations from a single, leptonically decaying W boson in 
the final state [23] . Here we simply note that the tv pro- 
duction is small compared to the main SM backgrounds, 
namely lU+jets and tt production, which have a similar 
signature. Thus we do not expect very stringent con- 
straints from the present single-top production searches. 

If X is a Majorana particle, then in priciple one can ex- 
pect a non- vanishing production of same-sign tops from 
the uu — > tt process; however, this amplitue is strongly 
suppressed by the smallness of m x . For m x = 2 GeV and 
mj > rrit the production of same sign tops is suppressed 



by (m x /m t Y 



10 with respect to the %-mediated 



stop anti-stop production (that is a small fraction of the 
SM tt cross section). Beside the (m x /m i ) 2 power sup- 
pression, the same-sign top cross section receives loga- 
rithimic enhencements both near threshold and at large 
dilepton invariant masses; however, these are compen- 
sated by the small probability distribution for the initial 
uu state. As a confirmation of this qualitative evaluation, 
we have perfomed a quantitative study of the same-sign 
top cross section for our preferred range of parameters 
finding cr(tt) T ovatron ~ 0.01 fb and cr(tt) [LH c 7TcV] ~ 3 fb: 
a level well below the present sensitivity. 



IV. FLAVOUR PHYSICS AND DARK MATTER 

A very significant constraint on Y q comes from the D°- 
D mixing amplitude, to which x° and t can contribute at 
one loop. Introducing the following low-energy effective 
dimension-six Hamiltonian 



(9) 



the one-loop induced contribution by our Lagrangian in 
(7) is 



1 



327r 2 m? 



(10) 



Using the bound \C* d \ < 1.2 x 10~ 3 TeV~ 2 [24], this 
implies 



2 At the LHC, the QCD production through gg initial state com- 
pletely dominates over the anomalous t-channel contributions for 
our choices of parameters. 



Yol% 



< 0.06 , 



(11) 



for our preferred values of mj and Y u . This is certainly 
a fine-tuned condition, although it is comparable to the 
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hierarchies exhibited by the Yukawa couplings within the 
SM. 

Another interesting consequence of our scenario are 
FCNC top quark decays to a light quark jet and miss- 
ing energy (i.e. t — > ux° X )- Neglecting the X ° mass 
dependence, the decay rate is given by 



r(t u X °x ) = 



\Y t Y u \ 2 m 5 t 
61447r 3 m4 ' 



(12) 



For our illustrative choice of parameters, the branching 
ratio might reach the level of 10 -3 , which could be within 
the projected LHC sensitivity [25]. 

Finally, it is interesting to note that stable fermions like 
X° with mass of a few GeV, annihilating to light quark 
pairs via effective dimension six operators have been pre- 
viously considered as dark-matter candidates [26] . In our 
case, assuming for the moment that x° is a Dirac fermion, 
the low-energy coupling of x° to light quarks induced by 
the t exchange leads to 



£Cff 



\Y U 



annih. 



4to? 



^WV(1-75)X°- (13) 



In this case the dominant contribution to the thermal 
annihilation rate of x° comes from the vector current 
part (x°7aiX°) 01 the above operator. Using the results 
of [26], and setting m~ t ~ 200 GeV and |F„| ~ 1.5 in 
(13) we find that the correct relic abundance of x° is re- 
produced for m x o ~ 2 GeV, that would perfectly fit with 
the collider data. Such contributions are however in some 
tension with existing direct dark-matter detection exper- 
iments [27]. Furthermore, a recent dark-matter search 
using the WMAP CMB [28] spectrum disfavors thermal 
relics of masses below 5 GeV. Therefore in order for X ° 
to have escaped these existing searches, its thermal an- 
nihilation cross-section has to be somewhat smaller than 
suggested by the dark-matter relic abundance measure- 
ments. This could happen if x° is °f Majorana type, 
so that its thermal annihilation cross-section is velocity 
suppressed. Alteratively, we could assume that x° re P~ 
resents only a fraction of the total relic abundance. A 
third possibility is to lower the x° mass below the sensi- 
tivity range of the present direct detection experiments, 
ignoring the cosmological bounds. If x° is a Majorana 
particle, then only the axial current (x°7At75X°) P ar t of 
(13) is allowed and current dark-matter searches (both 
direct and indirect) do not exclude a particle of this type 
with mass of a couple GeV or less (see e.g. [26, 27, 29]). In 
this case the correct dark-matter abundance can in prin- 
ciple still be obtained if %° is produced non-thermally 
as in the asymmetric dark-matter scenarios [30], for ex- 
ample through the out-of-equilibrium decays of heavier 
particles, not included in our low-energy model. 



pie analysed at the Tevatron is enriched by non-standard 
contributions coming from decays of a scalar top-partner, 
with the electric and color charge of the top quark. The 
top partners should be produced in pairs (iP), have a 
mass slightly above m t , and a large t — > tx° branching 
ratio, where x° is a SM singlet with mass of a few GeV 
or less (escaping detection) . With a proper tuning of the 
qtx° effective couplings, the additional ti pairs produced 
in this way could account for the large Ap B observed at 
CDF. 

The simple model we have proposed, where t and 
X° are the only relevant new light states, is fully con- 
sistent with present high-energy data and contains a 
dark-matter candidate. The model requires a non-trivial 
flavour structure, 3 such that t has large couplings to both 
X°tn a n d x° u Ri and a vanishing coupling to x° c B- How- 
ever, it is quite natural as far as the field content is con- 
cerned: a light scalar partner for the top is a natural 
candidate for a stabilization of the leading quadratic di- 
vergence on the SM Higgs mass term. In this paper we 
have pursued a bottom-up phenomenological approach 
and we have not analysed how this simple framework 
could be extended at high energies to be considered a vi- 
able extension of the SM; however, at this stage we see no 
obstacles to consider it as the low-energy side of a more 
ambitious supersymmetric model. 

Interestingly, this framework can soon be tested in 
more detail at the LHC. In particular: i) the large E™ ms 
of the sub-leading t — » ux° decay mode, and ii) the rise 
of the pp — > tP — > ti QCD cross-section at large m t f\, 
offer powerful tools to disproof or find evidences of this 
non-standard framework. 
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VII. ADDED NOTE 

After this work was completed, we become aware of a 
new anomaly reported by the CDF experiment, namely 
an excess in the production of jet pairs in association 



CONCLUSIONS 



Our knowledge of top-quark physics is still rather lim- 
ited. As we have shown, it could well be that the ti sam- 



s A flavour-violating structure of this type has also been postu- 
lated in [31], although in a different framework, and could be 
attributed to a specific flavour symmetry. 
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with a W boson [32]. In particular, CDF reports an ex- 
cess in the di-jet invariant mass (Mjj) distribution for 
120 < Mjj < 160 GeV. We note that our model pre- 
dict an excess over the SM prediction in the W + jj 
channel due to the sub-leading t — > u + x decay mode: 
tP -> ut (tu) + E miss -> ub (bu) + W-(W+) + E miss . 
We have simulated this decay chain and found that 
the resulting Mjj distribution is indeed peaked around 
140 — 150 GeV for our preferred range of parameters. 
However, the size of the cross-section is substantially 
smaller (about 1/3) with respect to the central value 



of the non-standard effect reported by CDF (after tak- 
ing into account the experimental cuts applied in [32]). 
Moreover, the Mjj distribution predicted in our frame- 
work is substantially broader, with a(Mjj) ~ 100 GeV 
with respect to the one reported in [32]. From these 
considerations we conclude that the inclusion of the tv 
production could partially improve the agreement with 
data also in the W + jj channel, although a quantitative 
evaluation of this improvement would require a more ac- 
curate evaluation of detector efficiencies and resolutions, 
especially on the jet variables, that is beyond the scope 
of the present work. 
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